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According to experimental data at GSI, the rates of the number of daughter ions, produced by 
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the nuclear K-shell electron capture (EC) decays of the H-like Pr 
ions, are modulated in time with periods Tec of the order of a few seconds, obeying an A-scaling 
Tec = j4/20s, where A is the mass number of the mother nuclei, and with amplitudes af c ~ 0.21. 
In turn, the positron decay mode of the H-like 142 Pm 60+ ions showed no time modulation of the 
decay rates. As has been shown in Phys. Rev. Lett. 103, 062502 (2009) and Phys. Rev. Lett. 
101, 182501 (2008), these data can be explained by the interference of two massive neutrino mass- 
eigenstates. In this letter we give a reply on the comments on our paper Phys. Rev. Lett. 103, 
062502 (2009) by A. Gal (arXiv: 0809.1213v2). PACS: 12.15.Ff, 13.15.+g, 23.40.Bw, 26.65. +t 



"GSI Oscillations" as interference of massive 
neutrino mass-eigenstates I2j 



■"^JTjie experimental data on the K-shell electron 
oaj>ture (EC) and positron (/3 + ) decay rates of the 
IJ^like heavy ions, carried out in the Experimen- 
L ral J Storage Ring (ESR) at GSI in Darmstadt Q- 
]4l] | place the following constraints on the theoreti- 
cal approaches for explaining of these phenomena: 
P")-the periods Tec of the time modulation of the 
04^-decay rates obey the A-scaling Tec = ^4/20 s, 
wiere A is the mass number of the mother ion, 2) 
the modulation amplitude for all observed decays 
' jC ~ 0.21 and 3) the /3 + -decay rates have no 
ic modulated terms. In addition all theoreti- 
cal approaches, describing the experimental data 
cjf|the EC and f3 + decay rates of the H-like heavy 
i^s] fl-Qj should explain the absence of the time 
[jK^dulation of the -EC-decay rates of bound atoms 
Wp m and 180 Re, measured in 0,0]. 

" "Afe has been shown in 0]-[ll[, the experimen- 
tal data 0-01, called the "GSI Oscillations", and 
the observations reported in 0, 0| can be explained 
only following the hypothesis of the interference of 
the massive neutrino mass-eigenstates \vj) 0,0] , 
defining the neutrino \v a ) with a lepton charge 
a = e, fj, or r as a coherent superposition of mas- 
sive neutrino mass-eigenstates \v a ) = Y^j^aji 
[12]. Neither the magnetic field of the ESR 
(see also fl3i] ) nor the mass-sp_litting of the H-like 



mother ions [9|, [ll| (see also [7[ ) enable to explain 
the experimental data The discussions of 

the "GSI oscillations" as the interference of mas- 
sive neutrino mass-eigenstates were proposed in 
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Recently [l5( , Gal has criticised the results ob- 
tained in [7| . The main aim of this letter is to clar- 



ify all problems, which have prevented Gal from 
accepting our results. 



Frequencies of "GSI oscillations" , caused by 
interferences of massive neutrino mass-eigenstates, 
should be inversely proportional to neutrino energies, 
i.e. to Q-values of EC -decays of H-like heavy ion s, 
but not to the mass number A of mother ions fldll 

The amplitude of the EC-decay m — > d + v e , 
caused by a Gamow-Teller 1 + — > + transition of 
the mother ion m from the ground hypcrfine state 
(ls) F= 1 into the daughter ion in the stable ground 
state and the electron neutrino v e , a coherent su- 
perposition \v e ) = J2j of massive neutrino 
mass-eigenstates Vj with masses mj, is a function 
of time t defined by 0] 

A(m -> dv e )(t) = ^2u ej A(m -> d,Vj)(t), (1) 
i 

where the amplitude A(m — > dvj)(t) of the m — ► 
d + Vj transition is calculated with the Hamilton 
operator of weak interactions given by 0, [1] 
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with standard notation 0, d, [l6| ■ It is equal to 
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where AEj = E d (qj) + Ej(kj) — M m is the energy 
difference of the final and initial state, M m is the 



mother ion mass, Ed(qj) and Ej{kj) are the ener- 
gies of the daughter ion and massive neutrino Vj 
with 3-momenta q } - and kj, respectively, A^gt is 
the nuclear matrix clement of the Gamow Teller 
transition m — > d and (ip\ s ) is the wave function 
of the bound electron in the H-like heavy ion m, 
averaged over the nuclear density [l6| ■ 

Suppose that in the -EC-decay m — > d + v e the 
momenta of the daughter ions qj, produced in the 
decay channels m — » d+Vj, are precisely measured. 
In case that the differences of momenta \q*i — qj\ 
are larger then the momentum resolutions | Sqj | , 
i.e. \qi — q*j\ 3> \Sq*j\, where <f, is a 3— momentum 
of a daughter ion in the decay channel m — > d + Vi 
for i ^ j, all decay channels m d + i/j are dis- 
tinguished experimentally and the EC-decay rate 
should never show time modulation. This agrees 
with the assertion pointed out in Ref . [ttT] . 

However, this is not the case with the "GSI oscil- 
lations" [Hi-El- The wave function of the detected 
daughter ion should be taken in the form of a wave 
packet, since the time differential detection of the 
daughter ions from the EC-decays with a time res- 
olution Td ~ 0.32 s introduces energy SEd ~ 2ir/Td 
and 3-momentum \6qa\ ~ 2ir/rdVd uncertainties, 
where Vd is the velocity of the daughter ion in 
the ESR 0. Such a smearing is described by the 
wave function &d(kj + q*j) 0. In the rest frame 
of the H-like mother ion the energy and momen- 
tum uncertainties are equal to 5 Ed ~ 2-K^/rd = 
1.85 x KT 14 eV and \Sq d \ ~ 2wyMd/ t^QeCi where 
7 = 1.432 is the Lorentz factor of the H-like 
mother ions 0, Vd = Qec/^cI is a velocity of 
the daughter ion and Qec is the C— value of the 
EC-decay, equal to the momentum of the daugh- 
ter, and Md is the mass of the daughter ion. For 
the EC-decay 14 °p r 58+ -» 140 Ce 58+ + v. the Q- 
value is equal to Q E c = 3348(6) keV 0. This 
gives \Sq d \ ~ 7.21 x lCT 10 eV. 

Due to energy and momentum conservation in 
every EC-decay channel m — ► d + Uj the energy 
and momentum of massive neutrino Uj are equal 
to 



!EC 



, \kj \ ~ Qec - 



,(4) 
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where Qec = — Afd- The differences of en- 
ergies and momenta of neutrino mass-eigenstates 
are 



Ei{ki) — Ej(kj) 
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In turn, w.y and fcy determine also the recoil en- 
ergy and 3-momentum differences of the daughter 



For two massive neutrino mass-eigenstates and 
the EC-decay of 140 P r 58+ we get 
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.40 x 10 _16 eV, 



3.27 x 10" n cV, 
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where we have set Am^ = 2.19 x 10 _4 eV 2 0. 
Since SEd 3> Wax an d |<5<7d| ^> |fc2x|, the daughter 
ions, produced in the two decay channels m — + d + 
i^! and m — * d + , are indistinguishable • This 
is the origin of the coherence in the EJC-decays 
m — > d + v e of the H-like heavy ions, measured in 
GSI 0. 

Thus, in GSI experiments the observed daugh- 
ter ion d is a nucleus with energy Ed{q) and 3- 
momentum q for all decay channels m — > d + i/j . 
The amplitude of the m —>■ d + v e decay reads 
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where AEj = E d (q) + Ej(fcj) - M ro . The EC- 
decay rate is related to the expression 



lim — - \ A (™ -» <^e)(i)| 2 = 3M m |M G x| 
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x^2E d (q)E j (k j ) $* d (ki + q) * d (Ay + q) 

x [2tt<J(A^) + 2irS(AE' j )] cos(wyt)}. (8) 

This expression reproduces Eq.® in our paper 
with the same frequencies = Am|/2M m . 

As has been mentioned in 0, the first term in 
Eq. © is the sum of the two diagonal terms of the 
transition probability into the states d + v\ and 
d + i/2, describing the incoherent contribution of 
massive neutrino mass-eigenstates, while the sec- 
ond term defines the interference of states Vi and 
Vj with i 7^ j causing the periodic time dependence 
with the frequency u>ij equal to 

ua = AEl - AEj = E d (q) + Ei(ki) - M m 

-E d {q) - Ej(kj) + M m = Ei(ki) - E 3 {k 3 ) = 
Am?, 

(9) 
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Thus, we argue that the interference term, pro- 
duced by a coherent contribution of massive neu- 
trino mass-eigenstates with a frequency inversely 
proportional to the mass number A of the mother 
ion, can be observed only due to energy and mo- 
mentum uncertainties, introduced by the time dif- 
ferential detection of the daughter ions ■ 

Due to the smallness of neutrino masses for the 
calculation of the -EC-decay rate we can take the 
masslcss limit everywhere except the modulated 
term U^Ugj cosfajt) 0. Since the 3-momenta of 
the massive neutrino mass-eigenstates and the 3- 
momcnta of the daughter ions differ only slightly 
from the Q-valuc of the EC decay, wc can set fcj ~ 
kj ~ k and \$ d (k + q)\ 2 = V(2ir) 3 5W(k + q), 
where V is a normalisation volume 0. Using the 
definition of the -EC-decay rate 
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we obtain the following expression for the time 
modulated EC-decay rate 



A-Ec(i) = A B c(l + a E c cos(w 2 ii)), 



(11) 



where u>2i = Am^ j2M m ^ a E c = sin26*i2 and Xec 
has been calculated in [16[. The .EC-decay rate 
Eq. (TTTj) is calculated for the matrix elements U e j of 
the mixing matrix U, taken at 6*13 = 00 (see also 
12j). In the laboratory frame the EC-decay rate 
is time modulated with a frequency ujec = W21/7. 
Thus, the period Tec of the time modulation is 
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For the experimental data on the periods of the 
time modulation 0-0 we get Amji = 2.19 x 
10" 4 eV 2 0. 



Coherence vs. incoherence in two-body electron 
capture - No interference terms in EC -decay rates of 
H-like heavy ions \Thl 

According to Gal the two-body K-shell 

electron capture m — > d decays, when only the 
daughter ions are observed, are driven by the com- 
plete set of orthogonal neutrino states \v a ) = 
with all lepton flavours a = e, fi and 
+ v n transition is 



T U 



a j I 

t. The amplitude of the m 
equal to [15[ 



A(m — > d v a )(t) 



3 



U a3 U ej A{m^ dvj)(t), (13) 



where the amplitude A(m — > dvj)(t) is defined by 
Eq.©. 

Since in GSI experiments neutrinos in the EC- 
decays of the H-like heavy ions are not detected, 
Gal proposes to define the probability of the m — > 
d transition as the incoherent sum of the squared 
absolute values of the amplitudes of the m — ► d+u a 
transitions [l5| 

P(m -+ d)(t) = \A{m dv a )(t)\ 2 = 

Ot 

= E E E U* al U* el U a] U ei A*{m -> d Vl ){t) 



a 1 j 

x A(m — > duj)(t), 
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where the index a runs over a = e, fi and r. Using 
the orthogonality relation for the matrix elements 
of the mixing matrix (l2^ 

E '*■'„'•.. .< ( 15 ) 

a 

one can arrive at the expression [l5| 



P(m -> d)(t) 
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|[/ e ,|V(m^d^)(i)| 2 ,(16) 



which contains no interference term. A similar ar- 
gument has been recently given by Yazaki [l8| . 

We want to point out here that the use of the 
complete set of neutrino wave functions \v a ) = 
J"V UXjWj) with all lepton flavours leaves room 
for the restoration of the interference terms in the 
rates of the m — > d transitions. 

The amplitude A(m — > d)(t) of the m — > d tran- 
sition we propose to define as a coherent superpo- 
sition of the amplitudes A(m — > dv a )(t) 

A(m -> d)(t) =^e- ,( »»A(m -> dv a )(t), (17) 

a. 

where we have introduced arbitrary phases tp a for 
neutrinos u a , which are responsible for the restora- 
tion of the interference term. The amplitudes 
A(m — > dv a )(t) are determined by Eq. (Ti"3"|) . The 
possibility to describe the amplitude A(m — > d)(t) 
by a coherent superposition Eq. (|17[) is obvious, 
since neutrinos v a are not detected. 

The rate of the m — ► d transition is related to 
the expression 
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where indices I and j denote the neutrino mass- 
eigenstates, the indices a and (3 run over all lepton 
flavours. The frequencies ujgj of the time modula- 
tion are defined by Eq.©. 

As has been remarked in 0] and discussed above, 
the first term in the r.h.s. of Eq. fTg)) corresponds 
to the decoherent contribution of massive neutrino 
mass-eigenstates, whereas the second one is caused 
by the coherent contribution of the decay channels 
m — ► d + v.j . 

The -EC-decay rates, measured in GSI experi- 
ments, take the form [lj 

Abc(£) = A B c(l + a E c cos(ujEct + 4>ec))- (19) 

In order to reproduce the correct value of the EC- 
decay constant Xec, calculated in [l6[ in the the- 
ory of weak interactions with massless neutrinos, 
we impose the following constraint on the phases 
of the neutrino wave functions 



i. 



(20) 



Setting the mixing angels #13 = and 6*23 = 7r/4 
and using the definition of the mixing matrix U 
12j , the condition Eq. (f2"0")) can be transcribed into 
the form 



5^ \J2U aj U ej e' 



(21) 



— i sin 2 26»i2 cos(w - (p Te ) + sm26*i 2 
2 V2 

X (cOSi^re — COS^e) COs26'i2 =1, 

where ip ae = tp a — tp e for a = p, r. 
The interference term is 

£ 2Re[j2u; i U* ei e i ^-^U od U ej e^ t 

= —= sin 26*12 (sin p^ e — snvp Te ) sin(u>2ii). (22) 
v2 

where we have used Eq. (f2"Tj) . 

Identifying the rate of the m-ni transition with 
the EC-decay rate Eq. (fTi?)) we get 
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sin 26*12 (sin tp^e — sin tp Te ) . (23) 



Thus, we predict that the phase and amplitude of 
the time modulated term of the EC-decay rate of 
the H like heavy ions should be universal and equal 
to 4>ec = — 7r/2 and a E c = ^75 811126*12 (sin ip^ - 
sin(^ Te ) (see Eq. ([23]) ) . respectively. 

For the experimental value of the mixing an- 
gle 6*12 = 34° [12fl and the modulation amplitude 
a ec = 0.21 [H-[3] we get the following system of 
equations 



sin (f^ e — sin p Te — 0.32 

1.75 COs(<^ e - tf Te ) + COSip^e - COS( y 9 re =0 



(24) 

defining the phase differences tp^ e and <p Te . The 
second equation we have obtained from the condi- 
tion Eq. (|2~lj) . One of the solutions of the system 
Eq. ([24|) is (fpe ~ 1.75 rad and tp Te ~ 0.73 rad. 

For the understanding of the physical origin of 
the phases tp a one can analyse, for example, the 
mixing matrix for neutrinos v a with definite lep- 
ton charges a = e, p and r as coherent superpo- 
sitions of massive neutrino mass-eigenstates. The 
wave functions of neutrinos u ai which we use for 
the calculation of the amplitudes of the EC-decay 
rates, can be written as \v a ) = J2j U^jWj), where 
the mixing matrix U* is defined by 
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(25) 

Here U* is the standard mixing matrix, calculated 
for the mixing angles 613 = and #23 = tt/4 flil ]. 
In such a representation of the mixing matrix, the 
phase differences tp^ e = tp^ — ip e and ip Te = tp T — tp e 
may have the meaning of Majorana phases, caused 
by CP-violation [H]. 

Since due to the condition Eq. ([2Tj) the limit of 
equal phases tp e = ip^ = tp T = (p, including tp = 0, 
of wave functions of neutrinos with lepton flavours 
a = e, p and r corresponds to the massless limit 
of massive neutrino mass-eigenstates rrij — > 0, be- 
ing equivalent to the vanishing of the mixing an- 
gle 6*12 — ► 0, the interference term vanishes for 
tp e = tpp = tp T = tp and, of course, for tp = 
as a partial case. 
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Final state interference in EC-decay rates is due to 
the neutrino magnetic moment [73 / 

The "GSI oscillations" cannot be induced by the 
magnetic moments of neutrinos, since in this case 
the experimental data on the EC and (3 + decays 
of the H-like ions should show a time modulation 
with equal periods. This contradicts 

Conclusive discussion 

We have shown that the appearance of the inter- 
ference terms in the EC-decay rates of the H-like 
heavy ions with periods, proportional to the mass 
of the H-like mother heavy ion Tec ~ M m but not 
the Q-value of the EC-decay, is due to overlap of 
massive neutrino mass-eigenstate energies and of 
the wave functions of the daughter ions in two- 
body decay channels m — > d + v\ and m — > d + , 
caused by the energy and momentum uncertain- 
ties, introduced by the time differential detection 
of the daughter ions in GSI experiments. 

We have shown that the idea that neutrinos with 
all lepton flavours contribute to the -EC-decays of 
the H-like heavy ions, which has been used by Gal 
to show the non-existence of interference terms, 
can be adopted for the derivation of the EC-decay 
rates with interference terms. 

We have pointed out that magnetic moments of 
massive neutrino mass-cigcnstates cannot be re- 
sponsible for the time modulation of the -EC-decay 
rates of the H-like have ions. The most important 
objective is that such a time modulation should 
be universal for EC and (3 + decay rates of all H- 
like heavy ions. This contradicts the experimental 
data . 
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